I. INTRODUCTION 



We obtain optical instabilities in all-optical bistable systems arising from competing co- 
operative pathways at low input light levels. In particular for three-level atomic systems 
in the A and V configuration interacting with two independent cavity modes, we identify 
the necessary conditions related to the incoherent pathways required to obtain instabilities. 
The instabilities arise when atomic states involved in the bistable transition are leaky and 
have substantial population, where the incoherent processes adversely affect the cooperative 
behavior of the atomic ensemble. 

Optical instabilities in atom-cavity coupled systems have been widely studied, particularly 



in the context of lasers and optical bistability 



iHsl . The phenomenon of all-optical bistability 



(AOB) is a cornerstone effect arising from cooperative atom-field dynamics. This phenomena 
has been extensively studied in order to realize all-optical switching, as well as nonlinear 
dynamical effects such as self-pulsing and chaos jo], Q]. We focus our attention to optical 
instabilities in AOB systems, particularly at low light levels. Traditionally the occurrence 
of optical instabilities arise in two well known regimes: one involves the atom coupling 
to multiple longitudinal modes of the cavity, and the other involves single cavity mode 
coupled to the atom and the instability occurs in the upper branch (also known as one- 
atom branch) of AOB jo-Sl, where cooperative effect is non-existent. The former class of 
instabilities, arising from large number of modes, result from an interplay of a variety of 
time scales leading to chaos and self-pulsing. The later effect involves intense intra-cavity 
field that saturates the collection of atoms in absorptive AOB along with optical pumping, 
thus does not involve any cooperative effects. In an effort to realize optical instabilities in 
AOB in the regime where cooperative effects play a central role, we propose coupling two 
adjacent transitions in a multi- level atom to bi-chromatic fields that experience independent 
feedback leading to competing cooperative behavior and the resulting optical instabilities. 
In our previous work, we have shown the existence of instabilities and negative hysteresis in 

nn 

three-level ladder (^) system with such double feedback AOB |9l. Il0l|. 

Earlier, three-level atom in the A-type configuration operated under Electromagnetic 
induced transparency (EIT) regime has been exploited to obtain instabilities using single 
cavity feedback [1V\. Recently controlled optical switching using double-cavity having a 



K-type multi- level system is proposed 12|. We seek to determine if instabilities can be 
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obtained universally in every few level atomic configuration. It appears that the nature of 
the incoherent process plays a central role apart from the cooperative effects. So far the 
three-level atom in the V and A-type configurations has not been explored in the context 
of double feedback AOB to observe optical instabilities. However, the three-level V-system 
AOB has been studied in the context of spontaneous emission induced quantum interference 
leading to variation in the AOB switching thresholds [isl as well as multistability fl^ . 
In absence of such quantum interference, the V system exhibits regular switching devoid 
of any optical instability 



15 



16| . In this paper we study both the A and V system in 



double-cavity AOB and are able to pin down the precise nature of the decay processes 
that could lead to controlled optical instabilities such as self-pulsing and chaotic output 



at low light leve 
communications 



s in the cooperative branch. Such optical chaos can be used for secure 
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as well as be exploited for synchronization of chaos 




FIG. 1: (a) & (b) Schematic of three- level atom in A & V configuration interacting with the 
two coherent fields having amplitudes Ei and E2 at frequencies ui and uj2, respectively. In case 
of V system, solid arrow from level |3) to |1) indicates incoherent pumping of rate r. (c) Two 
independent unidirectional ring cavities where the active medium contained within the length L 
interacts simultaneously with the two fields. 



The organization of the paper is as follows, in the next section we describe the necessary 
theoretical model foHowed by results and discussion associated with the A and V system, 
respectively. We conclude in the last section. 
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FIG. 2: (Color online) A system: (Left & right panel) Bistable response of output fields & x\ 
while the input field is varied and y\ is held constant. The red and green lines indicate stable 
steady states, while black line indicate unstable steady states which are physically inaccessible. 
The thick blue line indicates unstable steady states associated with instability. The parameters 
are |yi| = 0.5, Ci = 200, C2 = 2000, 71 = 1, 72 = 1, = 0.1, m = 1, K2 = 1, Ai = 0, A2 = 0, 
61 = 0,02 = 5. 

II. THEORETICAL MODEL 



The schematic in Fig. [T] indicates the three-level A and V configurations and the setup 
involving the double-cavity feedback, wherein two optical fields simultaneously interact with 
the atomic medium. We present the density matrix equations that govern the evolution of 
the atoms accompanied by the field equations that capture feedback in the slowly varying 



envelope approximation 



19| . We confine ourselves to the mean-field approximation, widely 



used in the context of AOB, thus ignoring any spatial dependence of the field within the 
cavity. The propagation of the two fields inside the active medium (atomic vapor) are chosen 
to be co-pro pag ating in order to eliminate the first order Doppler effect and the associated 
broadening [20 1. 

The two coherent fields having amplitude Ei and E2 at frequencies Ui and UJ2 couple to 
the atom in the A configuration along the transitions |2) o |1) and |2) ^ |3) respectively, 
are shown in Fig. [T]^a). The density matrix equations are obtained under the rotating- wave 



3 



approximation are given as 

dpii 
dt 

dpi2 



dt 
dt 

dp22 

dt 

dp23 

dt 

dp33 



271P22 - 2z/pii + iGip2i - iGlpu 
-(71 + l2 + i^ - 'iAi)pi2 - iGi{pu - P22) 

-iG*2pl3 

- [Z/ + i(Ai - A2)] pl3 + iGip23 - ^^2^12 (1) 

-2(71 + 72)p22 + iGlpu - iGip2i 

+iG2p32 - iG*2p23 

-(71 + 72 + «A2)p23 + iGlpis - iG2{p22 - P33) 

2upn + 272P22 + ^G2P23 - «G'2P32 



dt 

where Ai = CO21 — wi and A2 = ^23 — ^2 are the atomic detunings, Gi = ^21 ■ Ei/h, 
G2 = d23-E2/h are the Rabi frequencies, and 271 and 272 are the spontaneous emission decay 
rates from the level |2) — )■ |1) and |2) — )■ |3), respectively. The ground states decoherence 
arising from the non-radiative decay associated with the transition |1) ^ |3) is described by 
the decay rate 2z/. 

The three-level atoms in V configuration couple to two fields having amplitude Ei, E2 at 
frequencies Ui, U2 along the transitions |1) o |3) and |2) ^ |3), respectively, is shown in 
Fig. [l](b). The density matrix equations in rotating wave approximation are given as 
dpu 



dt 

dpi2 
dt 

dpi3 
dt 

dp22 

dt 

dp23 

dt 

dp33 

dt 



-271P11 - 2z/p22 + rp33 + iGip-ii - iGlpn 

-(71 +I2 + V + i(Ai - A2))pi2 + iGip32 
-iC^pis 

-{i' + 71 + ^A2)pi3 - iG2Pi2 - iGi{pii - P33) 

-272P22 + 2z/pii - zG';p23 + ^G'2P32 (2) 

-(72 + iA2)p23 - iG2{p22 - P33) - iGip2i 

271P11 + 272P22 - rp33 + iG'2P23 - iG2P32 

+iGlpr3 - iGipsi 



where, Ai = Ur3 — and A2 = ^23 — W2 are the atomic detunings, Gi = di^.Ei/h, 
G2 = d23-E2/h are the Rabi frequencies, and 271 and 272 are the spontaneous emission 



rates from the levels |1) — )■ |3) and |2) — )■ |3), respectively. Excited states decoherence 
resulting from a non-radiative decay pathways from |1) to |2) is indicated with a decay rate 
2v. An incoherent pump applied to transfer atomic population from level |3) to |1) plays an 
important role in the nonlinear dynamical behavior of V-system, whose strength is described 
with rate r 21 1. 




FIG. 3: (Color online) A system: (Top panel) The bifurcation diagram corresponding to the 
AOB curve shown in Fig. [2] (left panel) , where the stable steady states (red) and unstable steady 
states (blue) corresponds to normal switching and self-pulsed/chaotic output, respectively. The 
corresponding largest Lyapunov exponent is shown in the bottom panel. The parameters are same 
as in Fig. [2j 



Here, the optical fields having amplitudes Ei and E2 circulate within two independent 
ring cavities which provide sufficient feedback that leads to cooperative phenomena. The 
equations of motion of the fields consisting of information of boundary conditions as well as 
cavity feedback associated with the individual cavities are given as follows 



dxj 
dt 



Ki [-Xi{l + i9i)+yi + 2iCiPr, 



(3) 



where, the index 'i' refers to the fields 1 and 2 at frequencies ui and U2, respectively, and 
there involves normalized input fields i/i = dmn-El"^ /h'~^i\/T'i and the output fields Xj = 
dmn-EI"^^ /h'-)i^/Ti, where jg the amplitude of the input(output) field associated with 

the i^^ field that couples to the transition \m) ^ \n) whose dipole moment is dmn- The 
subscripts m and n denote the appropriate transition involving the atomic levels |1), |2) 
and 1 3), as shown in Fig. [1] (a) & (b). The other terms in the above equation are the cavity 
decay = cTj/£j, the scaled cavity detuning 6i = 5f/Ti and the cooperative parameter 
Ci = aiL/2Ti, where 6^ = (wf — uji)C/c is the normalized cavity detuning, L is the extent 
of the region which contains the active atoms, ai is the absorption coefficient associated 



with field at frequency coi, Ci is tlie total length of the corresponding cavity whose resonant 
frequency is u'^ and the transmission coefficient of the input /output mirrors is Tj. All the 
frequency units are normalized with respect to the spontaneous decay rate 72, unless specified 
otherwise. The derivation of above equations and the details of numerical techniques are 
along the lines discussed extensively in the Ref. {lo| . 
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FIG. 4: (Color online) A system: (Top) The temporal evolution indicating periodic self-pulsing 
and chaotic dynamics are shown in left and right panels, respectively. (Bottom) The associated 
frequency spectral density plots. The operating point \y2\ = 83.1 for periodic self-pulsing and 
\y2\ = 95.2 for chaotic behavior, and the other parameters are same as in Fig. [2l 

We have self consistently solved the above set of nonlinear atom and field equations in 
the steady state. We have undertaken a linear stability analysis, which allows us to map 
out the various regions of stability such as those involving stable switching, self-pulsing or 
chaotic output. The bifurcation diagram is also obtained in order to understand the various 
nonlinear dynamical features exhibited in these AOB double-cavity systems. 

III. RESULTS AND DISCUSSION 



The double-cavity AOB system interacting with three-level atoms exhibits a variety of 
nonlinear dynamical effects which seem to occur in all the three types of atomic configu- 
rations (A, V and S system) in appropriate regimes. It is such comparative analysis that 
provides insight into the optical instability and its dependence to decay channels. Earlier 
studies with three-level atoms have adopted a simplified approach wherein the intermediate 
level is adiabatically eliminated 22, l23|], and the effective interaction encompasses the two- 
photon transition without taking into the effects of the intermediate state. This assumption 



effectively rules out any optical instability 
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251 ]. Such simplified treatment leads to nor- 
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FIG. 5: (Color online) A system: Stability domain map indicating stable- fixed-point region (blue), 
the self-pulsing region (green), and the chaotic region (red) as the atomic detuning (A2) and the 
cavity detuning (^2) are varied. The other parameters are as given in Fig. [2] with \yi\ = 0.5 and 
\y2\ = 91.9. 

mal stable switching. In this paper, we focus our attention on nonlinear dynamical effects 
associated with the cooperative (lower) branch of AOB in the A and V atoms, which have 
not been presented earlier jo], [lo|. Our earlier study of S system in the double-cavity ex- 
hibits AOB, we found that the variation in the ratio of the cooperative parameter along the 
two adjacent transitions and the variation of the decays results in novel behavior such as 
self-pulsing, chaos as well as negative hysteresis. However, the A and V system are quite 
distinct in their behavior with regard to double-cavity AOB. 

It is important to bring out the difference between the instabilities presented here and the 
single mode instabilities that arise in the upper branch under strong field conditions, which 
are known to be well understood in the literature jsj-Sj. The conventional upper branch 
instability invariably results from an intense nonlinear atom-field interaction and a phase 
mismatched feedback arising from nonzero cavity detuning. The fast nonlinear response 
and slow optical pumping are known to be the two competing time-scales that results in 
self-pulsing. In contrast, our double-cavity instabilities arise in the lower cooperative branch 
of AOB without involving finite cavity detuning. The cooperative branch which intrinsically 
involves different physical processes is responsible for the instability. 

In case of the single feedback AOB system having a multi level atomic medium, the 
instability associated with the cooperative branch arises from leaky population of the lower 



atomic level involved in the AOB transition 26|]. The preparation of the atoms in the 



lower state such that they exhibit cooperative (or collective) phenomena is necessary to 
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observe normal switching. A leakage of atomic population from such a coherently prepared 
lower atomic level leads to instability. Such an effect is easily realized in the three-level 
ladder system with AOB along the upper transition. The leakage of population occurs via 
spontaneous emission into the ground state leading to the instability in the lower cooperative 
branch. Such an arrangement can also occur in the A system wherein incoherent decay 
coupling the pair of ground states is essential to obtain instability in the cooperative branch. 
This situation cannot arise in a V-system, as the AOB field directly couples to the ground 
state. Hence, V-system does not exhibit instability in the lower cooperative branch. 
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FIG. 6: (Color online) V system: (Left & right panel) Bistable response of output fields X2 &: xi 
while the input field 1/2 is varied and yi is held constant. The red and green lines indicate stable 
steady states, while black line indicate unstable steady states which are physically inaccessible. 
The thick blue line indicates unstable steady states associated with instability. The parameters 
are Ci = 200, C2 = 200, |yi| = 15, 71 = 1, 72 = 1, = 0.2, ki = 1,k2 = 1,Ai = 0,A2 = 4.5, 6I1 = 
0,02 = -2.5, r = 2. 

The various incoherent decay processes including spontaneous emission, relaxation of the 
ground state coherence as well as any incoherent pumping pathways to and from different 
atomic states play a critical role in dictating the nonlinear dynamical behavior of the double- 
cavity AOB system. Before taking the discussion further we recall some of the important 
results reported in our previous work jo], [lo|, which are pertinent in developing our under- 
standing of the current system. The three-level S system is coupled to a pair of coherent 
fields that simultaneously experience feedback in two independent ring cavities, apart from 
normal stable switching, the cavity fields exhibit either periodic self-pulsing or chaotic non- 
linear dynamical behavior in the cooperative branch. The system exhibits period doubling 
route to chaos. Another interesting effect is the presence of negative hysteresis at low light 
levels, whereas the conventional (positive) hysteresis continues to occur at higher incident 
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intensities. 




32 34 36 38 40 42 44 46 48 50 



FIG. 7: (Color online) V system: (Top panel) The bifurcation diagram corresponding to the lower 
branch of AOB curve shown in Fig. [6] (left panel), where the stable steady states (red) and unstable 
steady states (blue) corresponds to normal switching and self-pulsed/chaotic output, respectively. 
The corresponding largest Lyapunov exponent is shown in the bottom panel. The parameters are 
same as in Fig. [6l 



A. A - system 

The three-level A-system in a double-cavity configuration is presented here. We fix the 
strength of one of the input field say yi coupling the atomic transition |2) o |1) to a constant 
value, and then vary the other input field y2 (which couples the adjacent atomic transition 
|2) o |3)) to obtain the bistable response of y2 versus X2, as well as, y2 versus xi as shown in 
the Fig. O The variation of the cavity output field xi due to the change in y2 indicates the 
strong coupling and mutual dependence between the cavity fields Xi and X2 aided the atom. 
We have undertaken the linear stability analysis about the fixed points (steady state), and 
identified the unstable regions in the lower branch (indicated as thick blue line) as shown in 
Fig. O The usual unstable states associated with the negative slope of the S-shaped bistable 
response continue to exist (indicated as black line). The dynamical evolution of both the 
fields occur simultaneously and they mimic each other to a large extent. One can also 
obtain unstable states on the top branch by introducing finite cavity field detuning, akin to 
the instability obtained earlier 6|, |7|] . We further analyze these unstable states through the 
bifurcation diagram and these can be associated with either periodic self-pulsing or chaotic 
dynamics (see Fig. [2]&|1]). The stability of the periodic self-pulsing behavior is examined 
through the Floquet analysis whereas the existence of the chaos is confirmed through the 
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Lyapunov exponents (Ly) which take positive values corresponding to the chaotic regime 
(see bottom panel of Fig. |3]). As we vary the input field ?/2 the chaotic instability arises 
through a period doubling (PD) route as shown in the inset of Fig. [31 wherein at-least one 
of the Floquet multipliers (Fl) crosses out of the unit circle along the negative real axis j2^, 
specifically for I1/2I = 83.131 the dominant Fl is 1 and for \y2\ = 83.146 the dominant Fl is 
-1.0025. 




FIG. 8: (Color online) V system: The Limit cycle continuation from the Hopf point (H) at the very 
onset of instability on lower branch of AOB curve shown in Fig. [6] (left panel). The Limit point of 
cycles (LPC) as well as the period doubling point (PD) are indicated as cyan colored loops. The 
parameters are same as in Fig. [6l 

In order to obtain instability on the top branch, it is essential to have finite cavity 
detuning, we have chosen 62 = 5. In this regime the output amplitude of the chaotic 
oscillation extends from the lower branch to the top branch of the AOB response. The 
coexistence of chaotic and normal stable solutions occurs for higher input field strengths. 
The temporal evolution and the frequency spectrum of the cavity output fields for two 
input field strengths, as shown in Fig. IHfor two different |?/2| values, giving rise to periodic 
self-pulsing and chaotic behavior. The presence of the double feedback plays a crucial role 
in creating these instabilities. This instability is suppressed if the feedback for one of the 
fields is removed. The competition between atom-field nonlinear interaction along the two 
cooperative branches and the associated decay pathways within the atom play a critical 
role in obtaining these instabilities which are quite robust and occur over a wide parameter 
range. The parameter space associated with this system is extremely large with fourteen 
independent physical parameters involving decays, detunings and couplings. In order to steer 
the system systematically across these variety of instabilities we take recourse to stability 
domain maps, between two system parameters such as detunings 62 and A2 and are shown 
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in the Fig. |5l One can identify islands of chaos (red) within the periodic self-pulsing domains 
(green) as well as, regimes of normal switching (blue). 

B. V - system 

We would like to explicitly bring out the importance of the incoherent (decay) pathways 
within the atom that are crucial to the occurrence of instability. The absorptive AOB 
relies on collective absorption in the cooperative (lower) branch, in contrast to a saturation 
of independent atoms in the one-atom (upper) AOB branch. Thus, the preparation of 
the collective state in the atomic ground state is crucial for realizing regular switching in 
AOB. The incoherent processes that directly affect this collective preparation of the atoms 
leads to instability such as self-pulsing and chaos. A slow time scale can be associated 
with the decay of the initial preparation of the atom in ground state whereas nonlinear 
interaction of the atoms and the fields within the cavity provides the fast time scale, and 
a competition between these timescales leads to the self-pulsing behavior. We note that in 
the S and A configurations the intermediate state decay and the ground state decoherence, 
respectively, provide the incoherent decay pathways that turns the state responsible for 
collective absorption leaky, leading to instability. However, in the V system the shared 
ground state across the two transitions can not be leaky and does not lead to any instability 
in the lower (cooperative) branch of AOB. Essentially, it behaves like a pair of two-level atoms 
interacting independently with coherent cavity fields, coupled to a common ground state. 
Thus, in order to produce instability we again address the collective state and introduce an 
incoherent pump between the ground state to one of the excited states {i.e from |3) to |1) 
as shown in Fig. [T](b)). Introduction of this incoherent pathways leads to unstable steady 
states over the lower branch of AOB curve, and is shown in Fig. E] (indicated with thick blue 
line) and is essential for obtaining optical instability. The associated bifurcation diagram is 
given in the Fig. [TJ and again indicates the transition of the cavity output field from normal 
stable switching to chaos through a period doubling route, as shown in the inset. We also 
determined the Lyapunov exponents (see Fig. [7] bottom panel) to demarcate the range of 
input field \y2\ which lead to these instabilities. 

In general, with variation in the system parameter (I1/2I in the present case) the nonlinear 
dynamical self-pulsing behavior is expected to occur beyond the Hopf (H) bifurcation point 
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FIG. 9: (Color online) V system: (Top) The temporal evolution of periodic self-pulsing and chaotic 
dynamics are shown in left and right panels, respectively. (Bottom) The associated frequency 
spectral density plots. The operating point \y2\ = 38 for periodic self-pulsing and 1^2 1 = 47 for 
chaotic behavior, and the other parameters are same as in Fig. [6l 




FIG. 10: (Color online) V system: Stability domain map indicating stable-fixed-point region (blue), 
the self-pulsing region (green), and the chaotic region (red) as the atomic detuning (A2) and the 
cavity detuning (^2) are varied. The other parameters are as given in Fig. [U]with \yi\ = 15 and 
I2/2I = 50. 

that separates the normal switching from periodic self-pulsing. Interestingly periodic self- 
pulsing behavior coexists with stable switching for a range of I2/2I (see region around I1/2I = 
37 shown in Fig. [7]). To understand such coexistence we use the limit-cycle continuation 



from the Hopf point using the MATCONT continuation package [28|. The continuation 
diagram is shown in Fig. [8] and indicates that the limit-cycle continuation initially occurs 
in direction of decreasing |?/2| upto I1/2I = 35, and then turns around at a Limit-point 
of cycle (LPC), this is also known i.e Fold bifurcation point (indicated in cyan color). 
The dominant Floquet multipliers plotted on the complex plane indicate that periodic self- 
pulsing behavior (limit cycles) associated with continuation in the increasing and decreasing 
variation of \y2\ are stable and unstable, respectively. Therefore, all the three possible 
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Parameters 


Symbol 


Range 


Atomic detuning 


Al,2 


-10 ^ 10 


Cavity decays 


1^1,2 


o 1 


Cavity detuning 


^1,2 


-10 ^ 10 


Cooperative parameters 




100 o 1000 


Input field strengths 


|yi,2| 


o 150 


Non-radiative decay 


V 


o 1 


Incoherent pumping rate 


r 


o 1 



TABLE I: Ah the frequency parameters are scaled with 72 = 36MHz. 

behavior (normal switching and stable and unstable periodic self-pulsing) coexist in the 
lower AOB branch region within the interval |?/2| = [35:36]. Hence, depending on the initial 
condition the appropriate behavior will be observed. The period doubling route described 
in the above bifurcation diagram is also corroborated by the presence of the period doubling 
(PD) cycle indicated as cyan loop in Fig. [8l In the Fig. IH we present the temporal evolution 
of cavity output field corresponding to stable periodic self-pulsing and chaotic operating 
regimes, and the corresponding frequency spectrum in the bottom panels. The domain map 
plotted between the detunings 6*2 and A2 (see Fig. [10]) portrays the various regions such as 
normal switching (blue), periodic self-pulsing (green) and chaotic behavior (red) associated 
with cavity output. 

Before we conclude, we indicate the possible experimental system, where these effects can 
be realized, using Rubidium (*^Rb Di line) atomic vapor in the A and V-type configuration 
along the transitions 55'i/2(-F = 2) ^ 5Pi/2{F' = 2) o 5S'i/2(-F = 1) and 5P3/2 o 551/2 ^ 
5Pi/2, respectively. Considering temperatures of about 60°c one would obtain a number 
density of ~ 10^^ atoms/cm^ and with transmission coefficient T ^ 10^^ would result in 
cooperative parameter C ~ 1000. The input power levels can be varied from — 20 mW 
across a spot size of 100 fim to observe effects we have mentioned in this paper. The different 
values of cooperative parameter can be realized by changing either the number density of 
atoms or the transmission coefficient of the cavity. The range of parameters are explicitly 
enumerated in the Table [H The 4D3/2 ^ 5P3/2 ^ 4D5/2 transitions in ^^Rb offer promising 
implementation of these effects in the optical communication wavelength range of 1.5/xm 
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with degenerate fields 
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IV. CONCLUSION AND REMARKS 

We have demonstrated methods of generating periodic self-pulsing and chaotic instability 
using double-cavity AOB system having three-level atomic medium in the A and V configura- 
tions. These instabilities occur in the lower branch of the bistable curve and are intrinsically 
related to the competition of cooperative behavior of the atomic collection along the two 
adjacent transitions as well as the details of the incoherent pathways within the atom. The 
ground state decoherence is sufficient for the A system to exhibit instability, however, by 
introducing appropriate incoherent pathways (incoherent pump) we induce instability even 
in the V system double-cavity AOB. Hence, the incoherent pathway needs to deplete the 
collection of atoms involved in the cooperative effect. 
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